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I. Introduction 

Measured flow variables near the test section boundaries, used to guide adjustments of 
the walls in adaptive wind tunnels, can also be used to quantify the residual interference. 
Because of a finite number of wall control devices (jacks, plenum compartments), the 
finite test section length, and the approximation character of adaptation algorithms, the 
unconfined flow conditions are not expected to be precisely attained even in the ‘fully’ 
adapted stage [l],[2j. 

The procedures for the evaluation of residual wall interference are essentially the same 
as those used for assessing the corrections in conventional, non-adaptive wind tunnels. 
Depending upon the number of flow variables utilized, we speak of one- or two-variable 
methods [3]; in two dimensions also of Schwarz- or Cauchy-type methods [4]. 

The one- variable methods use the measured static pressure distribution at the test sec- 
tion boundary and supplement it with the far field representation of the model, estimated 
from its geometry and measured forces. 

The two-variable methods use measurements of static pressure and normal velocity at 
the test section boundary, but do not require any model representation. This is clearly of 
an advantage for adaptive wall test sections, which are often relatively small with respect 
to the test model, and for the variety of complex flows commonly encountered in wind 
tunnel testing. For test sections with flexible walls the normal component of velocity is 
given by the shape of the wall, adjusted for the displacement effect of its boundary layer. 
For ventilated test section walls it has to be measured by the Calspan Pipes, Laser Doppler 
Velocimetry, or other appropriate techniques. 

The interface discontinuity method , also described, is a ‘genuine’ residual interference 
assessment technique. It is specific to adaptive wall wind tunnels, where the computation 
results for the fictitious flow in the exterior of the test section are provided. 


II. Linear Flow Analysis 


Since the adaptive walls introduce only minor disturbances to the unconfined far field 
of the test model, the linearization of the potential equation near the walls is applicable 
as long as the flow remains subcritical there. 


(1) 


The governing equation for the disturbance velocity potential is 

,3V 3V 3V 
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where 0 = y/l — A/ 2 ^ and < 1 is the stream Mach number. 
The scaling of the streamwise coordinate, 



reduces Eq.(l) to Laplace’s equation, V 2 ^ = 0. 

The linear flow region where <f> satisfies Eq.(l) is shown schematically in Fig. la. It 
excludes the volume occupied by the test model, its viscous and transonic flow regions, and 
the wind tunnel exterior, where no real flow exists. The outer bounding surface, enclosing 
the test model, is expected to lie entirely within the linear flow region, off the viscous or 
nonisentropic flow at the walls. 

Using the principle of linear superposition, the disturbance velocity potential is split 
as [5] 

$ ~ "l~ > (3) 

where ^ TO is that due to the model in free air and <f> w is that due to wall interference. 

The model potential, <f> m > satisfies Eq.(l) in the infinite space outside the model and 
the adjacent nonlinear flow regions, Fig. lb. 

The wall interference potential, <f > w , is assumed to satisfy Eq.(l) in the entire test 
section interior, including the model and its nonlinear flow regions, as indicated in Fig.lc. 

This assignment of the singular and nonsingular parts as the effects of the model and 
the walls respectively is consistent with the concept of Green’s function for the Laplace 
operator. Accordingly, it is rigorous for an infinitesimal model, but only approximate for 
a finite-size model. 

The derivatives of <f> w are interpreted as disturbances to stream velocity components. 
They are usually evaluated at the model reference station or as averages over the model 
and interpreted as global corrections to stream Mach number |6] 

AA/„ = (L+X^Mi)M^“, (4) 

and to flow angles (in radians) 



From the spatial variations of these corrections over the model additional streamline cur- 
vature and buoyancy effects on model force data can be determined. 

In connection with adaptive wall wind tunnels, another type of the disturbance velocity 
potential is helpful: that corresponding to the ‘fictitious’ flow outside the wind tunnel. 
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The potential, denoted here by the symbol <j>, satisfies Eq.(l) in the exterior of the outer 
bounding surface, Fig.ld. The surface, separating the real wind tunnel flow and the 
computed exterior flow is also termed the interface. The aim of adaptation is to adjust the 
walls so that <f> and <f> constitute a single potential (f> m , continuous at the interface. There 
is a direct relationship between <p w and the difference <j> - <f> at the interface. 


A. One- Variable Method 

The method, due to Capelier, Chevallier and Dountol |7], is the most popular technique 
for the assessment of subsonic wall interference in wind tunnels with perforated walls. It 
retains the essential features of the classical wall interference approach [5], but replaces 
the idealized wind tunnel boundary conditions by 


( 6 ) 

dx 2 p ’ 

where C p is the measured boundary pressure coefficient. The control surface along which 
the pressure is measured should be some distance away from the wall, where the distur- 
bances of individual holes (perforations) are smeared out. The application of the method 
to test sections with slotted walls is more problematic as the flow becomes homogeneous at 
rather large distances from the walls, and the pressures measured directly on slat surfaces 
do not necessarily represent the averaged values. 

The axial component of wall interference velocity, 


ti fU 


d<j> w 

dx ’ 
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satisfying inside the test section 


/j2 d 2 u U) , d 2 u w d 2 u w 
fj 1 I ^ o — U ’ 


dx 2 dy 2 dz 7 

is obtained from its boundary values 

1 d<f> m 
Uw ~ ~2 p ~ dx 
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as a 


solution of the interior Dirichlet problem. The transverse velocity components, 


dcf> t 


d<h 


, w — - q — and w w — > 


can 


dy 

be obtained from u w by integrating the irrotational-flow conditions 


( 10 ) 


dv w du t 


and 


dw w du , 


dx dy dx dz 

along a path from the upstream end of the test section. 


( 11 ) 


The Dirichlet problem for Laplace’s equation is one of the best explored problems 
in mathematical physics and there are a large number of methods available to solve it 
numerically. A natural approach is to solve the problem in terms of the double layer 
potential |8], leading to a doublet panel method [9]. For simpler geometries, closed form 
solutions are obtainable using integral transforms [7] or the Fourier method [I0j-|l2j. 

The complex-variable treatment |7] of the two-dimensional problem leads, as pointed 
out in Ref. (4], to the Schwarz problem, consisting of determining an analytic function 
inside a domain from its defined real part on the boundary. Theory [13] shows that 
the integration of Cauchy-Riemann equations (irrotational-flow conditions) introduces an 
unknown imaginary constant, which needs to be specified in order to make the solution 
unique (specification of the upstream flow angle). 

The accuracy of the one-variable method depends greatly on the accuracy with which 
the free air potential <f> m can be predicted on the control surfaces [14], [15]. Since the 
far field of 4> m is normally evaluated using the measured model loading, subject to wall 
interference, the prediction tends to be more exact near a fully adapted stage. However, 
when compared to the relative size of the model, the adaptive test sections are usually 
much narrower than the conventional ones, so that the representation of flow near the 
walls in terms of the model far field may not be satisfactory. 

Another source of inaccuracy is the finite length of the test section and sparseness of 
the experimental pressure data. The boundary values of u w have to be interpolated or 
extrapolated over a complete boundary (closed or infinite), in order to make the Dirichlet 
problem soluble. The adaptive test sections, which are typically longer than the conven- 
tional ones, will have a slight advantage in this regard. 

The method can be used to monitor the reduction of wall interference corrections in 
the course of adaptation, but can also be incorporated into the adaptation algorithm [ 16]. 
Interference-free (unconfined) flow will be characterized by the vanishing boundary values 
of u„: 


Uu, = 0 on S. 


( 12 ) 


Compensation for errors of reference velocity or pressure [7], also called the autocor- 
rective property [15] or autoconvergence |17], is an important feature of the method. It 
applies within the limits of linearization and may be stated as follows: if the error of 
the (upstream) reference velocity Uoo is SUoo, then the perturbation velocities U - 
on the boundary will be offset by -SU 0 Q . The ensuing incremental correction, being of 
equal magnitude but opposite sign to the reference velocity error, will restore U as the 
reference velocity. 

For the one-variable method, working with measured boundary pressures p, the auto- 
corrective property can easily be verified by introducing the pressure coefficient 


C P = 


P - Poo 
\PooVl,' 
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and its error [6] 


P - (Poo + bp oo) 


6Cp kt-uX iH'iPooVi,) 

iU« 


C\ 


~ 2 


2 - (2 - Ml)C p 

SU oo 


U 0 


f/ D 


if 


c p ~o. 


From Eq.(9), considering d<t> m /dx invariant, the boundary value of u w is found to have a 
constant increment 


b U «ti ^ C | 


bU o 


(13) 


which is equal and opposite to the relative error of reference velocity. This incremental 
correction also applies also interior points since 


- , , Woo 

6u w (x,y,z) = - — — = constant 

U oo 

is a solution of Eq.(8) satisfying the boundary condition (13). There are no other possi- 
bilities, as the solution to a Dirichlet problem is unique. 

Besides compensating for genuine reference velocity errors, the autocorrective principle 
also establishes the correspondence between Uoq based on plenum pressure and actual 
stream velocity in ventilated test sections. 


B. Two- Variable Method 

Measurement of the static pressure and normal velocity distributions along the control 
surface opens the possibility of evaluating subsonic wall interference bypassing the model 
representation. The two-variable method is most easily applied to solid wall test sections 
where the walls can serve as control surfaces. 

Independent formulations of this concept using Green’s theorem are due to Ashxll and 
Weeks [18] and Cauchy’s integral formula (in 2D) due to J. Smith [ 4 ]. 

To describe the method, we introduce the position vectors of an interior point and a 
boundary point, 

ro = (zo,yo,zo) and r = (x',1/,2), (14) 


and denote by 


G(r 0 ,r) = - 


47r|r 0 - r 


(15) 


the fundamental solution (unit-strength source) for the three-dimensional Laplace opera- 
tor. 

Green’s second identity gives for a function <f> w harmonic in the test section interior 


<Mr 0 ) 



dG(r 0 ,r) 

dn 


G(r 0 ,r) 


dn 
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and for a function <p m harmonic in the test section exterior 


0 = 



d(7(r 0 ,r) 

dn 


G(r 0 ,r) 


d</>m(r) 

dn 


]dS. 


The differential and integral operations are taken with respect to the unsubscripted coor- 
dinates; 5 is the control surface (interface) enclosing the test section interior, and d/dn is 
the derivative along the outward normal to the control surface in the transformed space 

(•£ » y» ^)- 

Adding the above formulae and eliminating ip m from Eq.(3), we obtain the correction 
formula of Ashill and Weeks [18|: 


Mro) = fj t ,16) 

It expresses the interior value of the wall interference potential in terms of boundary values 
of the (total) disturbance velocity potential. 

Considering the entire space, Eq.(16) describes a sectionally harmonic function 4> w 
having a jump discontinuity <f> across the surface S. This differs from the more conventional 
representation of the wall interference potential by external singularities, where (f> w is 
continuous across 5 and harmonic everywhere except at the singular points. Of course, 
inside the test section both representations are equivalent. 

Physically, integral (16) can be interpreted as a surface distribution of doublets 


dG(r 0 ,r) 

dn 

and a surface distribution of sources 


with density 


4>(r) 


G(i*o,r) with density — ^ . 

on 

The normal component of disturbance velocity d<f>jdn can be measured directly, 
whereas the potential <f> has to be evaluated by a streamwise integration of the measured 
pressure coefficient, Eq.(6). 

Another possibility offers the integration by parts ! 19], converting the surface distri- 
bution of doublets into a surface distribution of (horseshoe) vortices 

n(r 0 ,r )~J dx' with density = ~^ C p( r )- ( 17 ) 


The far upstream and downstream terms are eliminated by the virtue of 


^( r o>r) -* 0 as x' — ► — oo and fl(i*o,r) — ► 0 as x' — ► oo. 

Taking in Eq.(16) the limit as i*o becomes a point of a smooth surface element, we 
obtain 


M r <>) = + //" kw 


aG(r »' r) - G (r„,r)^)] i S, 


dn 


dn 


To 6 S. 


(18) 
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A small circular neighbourhood of the singular point r 0 is to be taken out from the surface 
5 for the doublet integral; its contribution has already been accounted for by the isolated 
term l<£(r 0 ). There is no ambiguity concerning the source integral, as the contribution of 
a small circular element around the point 1*0 is zero. 

Another interesting relationship is obtained by substituting Eq.(18) in Eq.(3): 


</>m(ro) = \<t>{r 0 )-// s W r )^““ ~ G ( r °’ r ) Hrl d5, 


ro G S. 


(19) 


This formula, similar to that developed in Ref. (19], determines the boundary value of the 
free air potential, <f> m , from the measured boundary values of <}> and d4>/dn. Provided that 
the difference between the boundary values of <f> and <j> m is small, it may be possible to 
achieve <f> = <j> m in a single adjustment of the walls. Equation (19) will then play the role 
of a single-step convergence formula , a concept introduced in Ref.[20]. 

Alternative formulations of the correction method based on Green’s theorem are given 
in Refs.|2l] and (22), comparisons and accuracy aspects are discussed in Ref.|23|. Model 
representation, as shown above, is no longer required, but the sparseness of boundary data 
and incomplete test section boundary remain as a major source of inaccuracy. 

The specification of interference-free conditions in the two-variable method is straight- 
forward. Setting 4> w = 0 in Eq.(18) or <f> m = <t> »n Eq.(19), we obtain 


\Mro) = - f f s \<H*) 


aG(r °' r) C(r„,r)^>]<IS, r„6S, 


dn 


( 20 ) 


which interrelates the values of 4> and df/dn on the bounding surface of an adapted test 
section. 

The descent to two dimensions is accomplished by putting 


ro = (zo»!/o), r = (i # ,y), G(r 0 ,r) - — In |r 0 


and replacing the surface integrals by contour integrals. 

However, more readily applicable results are obtained using Cauchy’s integral formula. 
To illustrate this approach, we introduce the complex coordinate 


z = x' + iy = -z + iy 


( 21 ) 


and the complex disturbance velocity 


w 


dd> , d(f> . a 

(») = du(x,y) - tu(x, y) =0 1 r(x,y)- t^- (*,»). 


dz'~”' dy 

In accordance with Eq.(3), the complex disturbance velocity is decomposed as 


( 22 ) 


w(z) = U>m(z) + 


(23) 
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where w w is analytic in the test section interior and w m is analytic in the test section 
exterior. Applying the Cauchy integral formula to a counterclockwise oriented contour C, 
we obtain for an interior point zq 


w w (z 0 ) = 


'[ 

2m J c z - z 0 


and 


Adding the integrals 
mula [4j: 


o = -L / ^ dz . 

2*1 J c Z - Zq 

and eliminating w m from Eq.(23), we obtain Smith’s correction for- 
, v If w(z) 

= i7ij c 7^ iz ' ( 24 > 


expressing the wall interference velocity in terms of boundary value of the (total) distur- 
bance velocity. 

The Cauchy type integral along a curved path can be evaluated as indicated in Ap- 
pendix. Using Eq.(22), the components of the wall interference velocity are obtained as: 


Uw(*o,J/o) = ^Refu^eo)} and v w (x 0 ,y 0 ) = -Imltu^zo)}. (25) 


An example of wall deflections and wall pressures from the tests [24] of the 9-in chord 
CAST 10-2/DOA 2 airfoil in the 13-in by 13-in flexible-wall test section -of the 
Langley Transonic Cryogenic Wind Tunnel (TCT) is shown in Fig.2. The wall pressure 
distribution at the stream Mach number of 0.7 is subcritical as required. The downstream 
end of the integration contour was placed so as to cut off the last three pressure points, 
drifting away from the undisturbed flow conditions. The distribution of residual corrections 
along the wind tunnel axis, evaluated by the two- variable method, is shown in Fig. 3. The 
flow in the test section is not interference free, but considering the size of the model with 
respect to the test section, the corrections are certainly small. 

More detailed formulae, together with residual interference evaluated for the ON- 
ERA/CERT T2 flexible wall wind tunnel, can be found in Ref.[25). 

Considering the entire complex plane, Eq.(24) describes a sectionally analytic function 
w w having a jump discontinuity w across the contour C. This is obviously in contrast with 
the conventional representation of the complex interference velocity by external poles, 
allowing w w to be analytically continued across C, but only up to the location of the 
poles. 

The Cauchy-type integral (24) can be recast into the contour integral 


u>, 


>{zo) = [ 

Jc 


n{z) 


+ 


ct(z) 


ds, 


l c [2n{z 0 -z) 2jr(z 0 -.z)J 

where ds = \dz\ is the counterclockwise oriented contour length element. 


(26) 
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The integral can he interpreted as a line distribution of vortices 

J- with density ' 7 ( 2 ) = Re|«;(z) j = 9i( 2 ) ( 2 ^) 

2n(zo — 2 ) ill' 

and a line distribution of sources 

— _J with density o(z ) = -Im( w(z) -q~\ f = 

2*(z 0 -z) l \ dz \> 

where q t is the tangent component of disturbance velocity (positive in the counterclockwise 
direction) and q n is the normal component of disturbance velocity (positive in the direction 
of the outward normal). The correspondence with Green’s theorem approach is evident. 

The autocorrective property of Eq.(24) again applies (15) and is easy to verify. Starting 
with the reference velocity increment /if/*,, the boundary value of the ^-component of 
disturbance velocity 

U — U oo 


is found to have an increment 

V - ((/„ 

bxi = rr 


I 6Uoo) 


— u ~ — 


6U c 
Ur* 


From Eqs.(22) and (24) it follows for the increment of the complex disturbance velocity at 

an interior point 20 , 

. , ftSu [ dz c 

6ub4*o) = --T / - — — 

2 7T t J (j Z Zq 


Finally, from Eqs.(25) ^ 

^Uuj(^o» !/o) Su IT , 

U OO 

6v w {x 0 ,y 0 ) 0. 

A practical verification of the autocorrective property is shown in Fig.4. The reference 
Mach number of our example in Fig.2 was tentatively changed from 0.700 to 0 695 and 
the wall pressure coefficients, used as input for the residual interference calculation, were 
recalculated accordingly. Comparing Fig.4 with Fig. 3, we note that the resultant Mach 
number correction curve is displaced by 0.005 in the positive direction, so that the corrected 
Mach number is again the same. The angle of attack correction, as expected, is not greatly 
affected by the change of the reference Mach number. 

Correction formula (24) is closely related to wall adaptation criteria for two- 
dimensional testing. In the limiting process, as Zo becomes a point on a smooth segment 
of the contour C we obtain 


w, 


■ I 2 ”) = -2 “<*0 ) + ^ / c 7TT„ dz ’ 20 6 C ' 


(29) 
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where the (singular) integral is to be interpreted as Cauchy’s principal value. 
Substituting Eq.(29) in (23), we find 


rv. 


. x 1 , . If tv(z) . 

>( zo ) = zW (Zo) - / - L ~ L - dz , z 0 <E C, 

2 27Tt Jc z- z 0 


(30) 


which is the limiting case of the formula given in Ref.[26|. It determines the boundary 
value u/ m of the complex disturbance velocity due to the model in free air, in terms of 
the measured values w. This result proves again that the model representation in the 
two-variable method is, in theory, superfluous. However, for incomplete boundary data an 
independently estimated far field of w m may conveniently be used to aid the interpolations 
and extrapolations. 

Equation (30) may also be used as the two-dimensional single-step convergence for- 
mula ; the case of straight line boundaries can be found in Refs.|20) and [27]. 

Setting w w = 0 in Eq.(29) or w m = w in Eq.(30), we obtain the interference-free 
condition 

1 , . If u>(z) 

2"' ( * o) = -2 7ij c 7±Z i2 ' *> SC < 31 > 


in terms of the complex disturbance velocity on the boundary. The factor | was left 
uncancelled, to emphasize the connection with the three-dimensional condition, Eq.(20). 

Considering straight line boundaries at y = ±|, we obtain in terms of disturbance 
velocity components 


, h. 

1 

r viz ' ± ' ) <“ 

u *o,±-) = 


2 


J — oo ^ *0 

v(i 0 ,±J) = 

± f ~. 

f 

2 

n J 

La, x ~ x o 


(32a) 

(326) 


These ‘compressible-flow’ versions of Hilbert’s transforms, introduced by Sears (lj as func- 
tional relationships between two velocity components, define unconfined flow in a two- 
dimensional test section. 


C. Interface Discontinuity Method 

This residual interference method, closely related to the two-variable method, utilizes 
exterior flow calculations. The general idea, as proposed by Sears and Erickson [28) is 
essentially this: the flow field is considered to consist of an experimental inner region 
joined at an interface to a computed outer region. If the computed outer flow satisfies 
the unconfined flow conditions and matches along the interface the inner flow, then the 
combined flow field is continuous, representing unconfined flow around the model. The 
matching error, or discontinuity, provides a measure of the residual interference. It can be 
quantified by removing the discontinuity by a surface distribution of singularities. These 
singularities do not disturb the unconfined flow condition in the outer region, but do 
introduce velocity perturbations at the position of the test model, which then can be 
interpreted as the usual wall interference corrections. 


I 
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As for the two-component method, Green’s theorem will give us a quick answer as to 
what the suitable singularities and their densities should be. Selecting 1*0 to be an interior 
point, we obtain for the function </>, representing the disturbance velocity potential of the 
fictitious flow in the exterior region 


0 = /// (r) 


dG(r 0 ,r) 

dn 


C ( r o>r)-j n -\dS. 


Subtracting it from Eq.(l6), we obtain the interior value of the wall interference potential 
in terms of the differences of the interior and exterior flow potentials and their normal 
derivatives along the interface: 


Physically, integral (33) can be interpreted as a surface distribution of doublets 


^! r ^ r ! with density 


[<t>{ r) ~ </>(r)] 


and a surface distribution of sources 

G(r 0 ,r) with density 


, d<f>(r) _ d<ft(r) , 
^ dn dn ' 


The potential <f> is obtained by solving an exterior flow problem (CFD), but <p w is 
obtained by a surface integration, as in the two-variable method. 

The exterior flow can be calculated as a solution of a Neumann problem, satisfying 
the boundary condition 

= reS , (34) 

dn dn 

where d<J>(r)/dn is the normal component of disturbance velocity on the interface. Integral 
(33) then reduces to the distribution of doublets, 


<M r o) 


ffim - <£(*•)] 


dG(r 0 ,r) 

dn 


dS. 


(35) 


Alternatively, the exterior flow can be calculated as a solution of a Dirichlet problem, 
satisfying the boundary condition 


J(r) = (£(r), r€ 5, 


(36) 


integral (33) reduces to the distribution of sources, 


<Mr 0 ) 



d?(r) 

dn 


]G(r 0 ,r)dS. 


(37) 
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The latter approach has recently been described by Rebstock and Lee |29|. 

Finally, if the walls are adjusted to satisfy the conditions (34) and (36) simultaneously 
(a perfect match), then from Eq.(33) 


<M r o) = 0, 


indicating that the flow inside the test section is interference free. The conditions of (low 
tangency and equal pressures along the interface imply that the desired interface is a 
stream tube. This streamlining principle for an adaptive wall test section, introduced by 
Goodyer |2|, is of course quite general and not just restricted to linear subsonic flow. 

The Cauchy integral approach, applicable to two-dimensional flow, proceeds along 
the similar lines. Considering the complex disturbance velocity ui of the fictitious flow, 

analytic in the exterior region and vanishing at infinity, then for an interior point z n it 
follows 



dz. 


Subtracting it from Eq.(24), we obtain 



(38) 


This Cauchy type integral can again be evaluated as described in Appendix. 

If the norma] component of disturbance velocity is continuous across the interface, 


&»(*) = 0«(*)> *ec. 


then from Eqs.(27)-(28) 


(39) 


l M*o) ~ ~ 9»(*)] 2*(z 0 Z z ) tls ‘ (40) 

The wall interference velocity is represented by contour distribution vortices, whose 
density is equal to the discontinuity of the tangential component of velocity. 

Conversely, if the tangential component of disturbance velocity is continuous, 

&(*) = ?<(*), zeC, ( 41 ) 

then 

M*>) = j -{(<,„ M - »„(*)] ^1—^ ds. (42) 

The wall interference velocity is represented by contour distribution sources, whose den- 
sity is equal and opposite to the discontinuity of the normal component of velocity. 

The single-step convergence formula of Judd, Wolf, and Guodyer [30| can be derived 
from Eq.(38) by taking the limit as as z 0 becomes a point of interface C, by analogy with 
Eqs.(24) and (29), and eliminating w w from Eq.(23). 
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Appendix 

The Cauchy type integral, Eq.(24) or (38), is easily evaluated by using a technique 
from Ref. (31). Approximating the contour C by line segments, the integral 



(43) 


reduces to the sum 

w w {z 0 ) = AjU>u,(zo), (44) 

} 

where 

i r z >+ i fiz) 

a '“"( 2 o) = i„i 4 2 ” 0 * < 45 ) 

is the contribution of the j-th segment. 

Assuming a linear variation of the density function / between the segment end points 
Zj and 2y + j: 


m = /, + J’ {*-*>) 

z j + 1 z j 

— [il l - h ( 2 - Z q) + /j4 

2 >+ 1 2 j 

and substituting it in Eq.(45), we find 


2o - 2> 


z j + i 


l- - /> 


2p - 2>+l 
2 J + 1 — Z J 


A ' u ’" (zo) = + Si 


/>+. * i" - 1, 

Z ) + 1 


Zp - 

2 rf i “ e j J 


In 


*>+] 


Zo 


Z| Zq 


(46) 


(47) 
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Fig- 1 Linearized flow regions. 
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Fig. 2 Wall deflections and wall pressure coefficients; 9-in chord CAST 
10 airfoil in the 13-in x 13-in test section of NASA TCT, Moo = 0.700, a = 1.20°, 

C N = 0.50. 






